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Abstract. Electromagnetic field fluctuations are responsible for the destruction of 
electron coherence (dephasing) in solids and in vacuum electron beam interference. 
The vacuum fluctuations are modified by conductors and dielectrics, as in the Casimir 
effect, and hence, bodies in the vicinity of the beams can influence the beam coherence. 
We calculate the quenching of interference of two beams moving in vacuum parallel to 
a thick plate with permittivity e(cij) — eg + i'^T^cr /uj. In case of an ideal conductor or 
dielectric (|e| = oo) the dephasing is suppressed when the beams are close to the surface 
of the plate, because the random tangential electric field i?t, responsible for dephasing, 
is zero at the surface. The situation is changed dramatically when eq or a-re finite. In 
this case there exists a layer near the surface, where the fluctuations of Et are strong 
due to evanescent near fields. The thickness of this near - field layer is of the order of 
the wavelength in the dielectric or the skin depth in the conductor, corresponding to a 
frequency which is the inverse electron time of flight from the emitter to the detector. 
When the beams are within this layer their dephasing is enhanced and for slow enough 
electrons can be even stronger than far from the surface. 
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1. Introduction 



Quantum electromagnetic (EM) field fluctuations are well known as being responsible 
for the Casimir forces, see for example jT]. Less known is the role of these fluctuations 
in destructing electron coherence. In weak localization phenomena in solids EM 
fluctuations are one of the dephasing mechanisms of conduction electrons [2], see also 
|3j. The interference of vacuum electron beams, observed experimentally [HE], is also 
quenched by EM fluctuations [HI, see also [7j. These two papers consider EM fluctuations 
in vacuum or when ideal conductors are present in the vicinity of the beams. The role of 
dissipation was discussed in Ref . jH] . Based on physical arguments, the decoherence was 
related to the deceleration of an electron from the beam due to the Ohmic dissipation 
of the current produced in the metal by the image charge. 
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The aim of this paper is to extend the calculations of Refs. jSl El to the case 
where the beams are close to dissipative bodies and to consider in detail the experiment 
geometry when two interfering beams move in vacuum parallel to a thick infinite plate 
with permittivity e^uj) = + i Ana/u. Calculations of the dephasing factor in this 
geometry demonstrate the crucial role of dissipation. If the plate is an ideal conductor, 
a = oo, the fluctuations of the tangential electric field Et, which are responsible 
for beam dephasing in this geometry, are suppressed near the plate surface because 
of the boundary condition Et = at the surface. However, when a is finite, very 
strong fluctuations of Ef exist near the plate surface, within a layer of the order of 
the skin depth. These near-field fluctuations dramatically enhance the beam dephasing. 
Unexpectedly, a similar effect exists also near a lossless dielectric with high permittivity, 
(T = 0, eo ^ 1, within a layer of the order of the wave length in the dielectric. 

The paper is organized as follows. In SecElwe present the dephasing factor in 
terms of the EM field correlator D in the case of no dissipation, Eq.©, and give reasons 
why the quantum Langevin equation for the EM field has to be used when dissipation 
is present. In Secj^Jwe derive the Langevin equation and prove that the expression of 
K in terms of V is valid in the case of dissipation too. In this section we present also 
the relation between V and the EM field retarded Green function Q, which is used to 
calculate V. The above mentioned special geometry is considered in Sec|31 where K is 
given as an integral, Eq. ()H2|l . over wave vectors and frequencies, containing the spectral 
density of the EM field fluctuations {EtEt)]^^;, and the spectral density |(l)k(^P of the EM 
field radiated by the beam electrons. {EtEt)kuj is calculated in SecEland SeclHl where 
Eqs. fj53|) and ()54|) demonstrate the enhancement of fluctuations due to near fields. In 
SecEl we present a model for |(l)ka;P and calculate explicitly as a function of the 
distance of the beams from the plate d and the electron velocity v (see Eqs. ()6H) . (jU^ 
and ()63|) and the text which follows). It turns out that the dephasing enhancement 
due to near fields is appreciable when v and a are not very large. In Sec|H]we discuss 
the relation between beam dephasing and beam EM radiation. The Appendix contains 
some calculation details. 

2. Beam dephasing 

If one ignores the interaction of the beam electrons with the EM field, the number of 
electrons measured in the interference experiment is n = |'?/;ip + |'?/'2p + 2Re('?/'iV'2)) where 
ipi and 1^2 are the wave-functions corresponding to the coherent motion of the electrons 
in beams 1 and 2, and n is calculated at the detector position. The interaction with 
the EM field does not affect the squares and \ip2\'^ (since it does not change the 
number of electrons in the beams), but the product ipiip2y responsible for the interference 
pattern, is multiplied by a factor e*'^e~^ with real and positive K. The first factor 
only shifts the interference pattern in space, while the second one reduces the amplitude 
of the interference oscillations (compared to the background I'^/'iP + |^/'2p) and describes 
dephasing. 
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To calculate the strength of dephasing we use the "trace of the environment" picture 
[2]. At t = to, when the electron is emitted from the source, the environment is in state 
\to)- While moving, the electron interacts with the environment and perturbs its state. 
When the electron moving in beam 1 arrives the detector at time ti, the environment 
evolves due to this interaction to state \ti). In a similar way one defines the state 1^2)- 



According to the "trace of the environment" picture e*''e 



-K 



{t2\tl). 



One can present the final states of the environment in terms of evolution operators, 



\ti) = uM, 



Ui = T exp 



dtHiit) 



to 



(1) 



where T means time ordering and Hi{t) is the interaction of the electron in beam 1 
with the environment. Hi{t) is in the interaction representation, i.e. sandwiched with 
evolution exponents exp [— iTit/^] containing the beam electron Hamiltonian and the 
environment Hamiltonian. In a similar way one defines U2 in terms of H2{t) and finds 
(^2|^i) = (^o|^2~^^il^o)- For an EM environment, choosing a gauge with zero scalar 
potential, we have 

Hi{t) = --^ j rfrji(r,t)A(r,t), (2) 

where ji electron current density operator for the electron in beam 1 (sandwiched 
with the evolution exponents containing the beam electron Hamiltonian) and A the 
vector potential operator (sandwiched with the evolution exponents containing the EM 
environment Hamiltonian). H2{t) is defined similarly with the current j2. 

In this approach one assumes that at the initial moment to the electron source and 
the environment are un-correlated. It is also assumed that the renormalization of the 
bare electron parameters due to the interaction with the EM environment [7] does not 
influence substantially the dephasing phenomena. 

To proceed we assume, following [H] , that the current is a classical quantity. When 
there is no dissipation in the EM environment, its Hamiltonian is simply the EM field 
Hamiltonian and the EM field can be quantized expanding it in normal modes. It is 
well known that in this case the commutator [A(r, t), A(r', t')] is an imaginary c-number, 
and due to the classical nature of the currents ji and j2 the commutators of Hi(t) and 
H2{t') have the same property. Because of this property the time ordering affects only 
the phase of the evolution operators P, and one can obtain 

i 

h 

if Hi is defined to be zero for t < and for t > ti. The phase 
commutator [Hi(t), Hi{t')]. Defining H2{t) in a similar way, we have 

i 

h 

where the additional phase x contains the commutator [ifi(t), 7^2(^0] • Averaging this 
over I to) we find 

i 

he 



Ui = e''*"Vi, Vi = exp 



dtHi{t) 



(3) 

1 contains the 



-<^2)gixexp 



dt{Hi{t) - H2{t)) 



(4) 



(^21^1) 



e^e 



-K 



-4>2+x) 



exp 



dt / cir ji2(r,t)A(r,t) 



(5) 
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where ji2(r, t) = ji(r,t) — j2(r, t) and (...) means average over \to). When the initial 
state of the environment is an equihbrium state with temperature T, the average means 
a thermal average {■■■)t- 

The second important property of A(r, t) in the case of no dissipation is that it 
is a Gaussian operator with respect to thermal averaging {■■■)t- After expanding A in 
normal modes this property follows from the relation P 

(exp(a*a^ — aa) )t = exp 

where is the bosonic operator creating a photon in some normal mode, n = {a^a)T is 
the occupation number of this mode, and a is a complex number. Using the Gaussian 
properties of A one can perform the thermal averaging in Eq.(jSI) and obtain 

K = I dtdt' I drdr'jUr,t)f,,{r',t'){A^{r,t)Ap{r',t'))T, (7) 

where a, (3 = x,y, z. If one defines the thermal correlator 

V^p{r,r';t-t') = ^{A^{r,t)Ap{r' ,t') + Ap{r' ,t')A^{r,t))T, (8) 
the final result is 

K = j dtdt' j rfrdr'j- (r, t)jUr', t')V^p{j, r'; t-t'). (9) 

It was obtained for T = in Ref . |6| and for T 7^ in Ref . [7] . We derived it in a different 
way to emphasize the two assumptions under which this result is valid (for classical 
currents), namely: (i) the commutator of the field operator A(r,t) is an imaginary 
c-number and (ii) A(r, t) is a Gaussian quantity with respect to thermal averaging. 

When dissipation is present, the EM environment Hamiltonian includes not only 
the EM field, but also the electrons in the absorbing bodies and their interaction with 
the EM field. If the field operator A is defined as sandwiched by evolution exponents 
containing the EM field Hamiltonian only, it has to be considered as a random quantity 
due to the influence of the dissipative electron system in the absorbing bodies. These 
electrons are the thermal bath, whose temperature defines the temperature of the EM 
field. Being a random operator, A obeys the quantum Langevin equation, where the 
effect of the dissipative electrons is simulated by a random force. We will show in what 
follows that the crucial properties of A(r, t) used to derive Eq.ijH} are valid also for the 
random vector potential operator, and hence Eq.ijH} is valid when dissipative bodies are 
present. Note, that in case of dissipation normal modes of the EM field do not exist, 
the EM field can not be quantized in the usual way, and this is why one is forced to use 
the Langevin equation approach. 

3. Quantum Langevin equation for the EM field 

A quantum Langevin equation for the coordinate operator g of a particle moving in 
potential V{q), derived in Ref. ^J, can be written in terms of the particle Lagrangian 



-{{a*a^ — aaY)T 



exp 



— \a\ 



1 

n -\ — 
2 



(6) 
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L = + V{q) as 



+ r rft'7(t-t')^(o =i'w- (10) 



dt dq{t) dq{t) 

The kernel 7 is responsible for the "friction" produced by the environment, which is 
a thermal bath at temperature T, and the operator F{t) is the random force. The 
statistical and commutation properties of the random force are defined by the dissipation 
kernel 7. Namely, F{t) is a Gaussian stationary random process with {F) = and a 
correlator 

l{F{t)F{t') + F{t')F{t)) = ^ r doo expl-tooit -t')]nwcoth^Re^{uj), (11) 
while the commutator of the random force is an imaginary c-number. 



2 f°° 

[F{t),F{t')] = - / dujexp[-iuj{t-t')] hiuRe-f{iu). 

i-00 



(12) 



(Note, that this Langevin equation is equivalent to the well-known approaches used by 
Feynman and Vernon ^1] and Caldeira and Legget [T2]). 

Consider the classical Maxwell equations for the long- wave EM field, rotE = — H/c 
and rotH = D/c + (47r/c)j, where D is the displacement given by 

D(r, t) = E(r, t) + [ dt' x(r, t - t')E{r, t'), (13) 

and j is the external current density. Note that fields entering the above macroscopic 
equations are averaged over a volume AV^ = (AL)^, where AL is large compared to all 
relevant microscopic lengths, but small compared to the wave-length of the EM field. 
With the gauge E = — A/c and H = rot A the first Maxwell equation is satisfied and 
the second gives an equation for the vector potential 

1 •■ IT- A-K 

^A(t) + rotrotA(t) + — / dt'A{t')xit - t') = — j(t). (14) 

J c 

Starting from the EM field Lagrangian one can prove that this equation can be 
considered as the quantum Langevin equation for the random field operator A, if j 
is the appropriate random force created by the thermal bath of dissipative electrons. 
[jaAV/c plays the role of F and x^^/47rc^ plays the role of 7.) The correlator of this 
force is known from the fluctuation-dissipation theorem [TH] , 

^ {Ur,t)jp{r',t')+jpiv',t')Ur,t)) = j du e-^-(*-*')(j«(r)j^(r')). (15) 

with 

{3a{r)jp{r'))^ = 6ai36{r -r')^ ^^coth^ Ime(r,aj). (16) 

Comparing this correlator with Eq.()ll|). we can find from Eq.()12|) the commutator of 
the random currents, which turns out to be an imaginary c-number, 

[j4r,t),jp{r',t')] = 6o.p6{r-r')^^ j du;e— (*-*V Ime(r, u;). (17) 
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(The quantum Langevin equation for the EM field was considered also in Ref.^l], but 
in a form not suitable for our problem.) 

The retarded Green function, corresponding to Eq. lfT^ . obeys 

^Gaxir, r'; t) + rotrot„^e?;3A(r, r'; t) + ^ j dt'Q^xir, r'; t')x{t - t') (18) 
= -47r5„A5(r-r')5(t), 

with the condition ^(t) = for t < 0. (The rot operators are defined in terms of the 
antisymmetric tensor 5aap as rot^/j = ^acr/jVo- and rotroto-/? = VqV/j — dapV^-) Using 
this Green function one can calculate the random field operator 



A„(r,t) = --^j dt' j dr'ga,3ir,r';t-t')j^{r',t'). (19) 



Two important consequences follows from this relation. First, as the commutator of the 
currents is an imaginary c-number and the Q is real, the commutator of the random 
field operators is also an imaginary c-number. Second, as the current is a stationary 
Gaussian process and Q depends on t — t', same is the random field operator. Since 
these two properties of the field operator were crucial for deriving the dephasing factor 
as given by Eq.® for the case of no dissipation, we proved thereby that this result is 
valid also when dissipation is present. 

The last equation allows also the calculation of the correlator and the commutator 
of the field operator. The correlator is known ^3] to be related to the Green function 
defined by Eq.dlHI), 



with 



V^p{r,t;r\t') = j duj e"^-(*-*')(A„(r)A^(r')). (20) 

{A^{y)Ap{y'))^ = -- coth ^ lmg^(,{iu; r, r'), (21) 

where the Green function in the frequency domain is defined as 

Qpx{v, r'; t)= j ^e-^-%,{u;; r, r'). (22) 

The definition of g{uj) corresponds to that in [16j. The Green function is symmetric, 
Qi3x{u!;r,r') = ^^/^(t^; r', r), and it follows from the properties of e{r,u!) that 
Ql3\{uj] r, r')* = Q/3\{—u!] r, r'). One can prove the following important integral relation 

— / rfroIme(ro,t^)^Aa(t^;r,ro)^A'a(t^;r',ro)* = -47rIm^AA'(^^;r,r'), (23) 

which is an obvious generalization to 3D of the ID relation given in ^3]. (It can be 
proven using the following Green theorem: f dr (p rotrota^-?/' = f dr Toto-a<P rotg-/?'?/'-)- 
Using this integral relation one can obtain Eq. ()21|) and also calculate the commutator 

2h 

[A,(r)^, Af,{rX,] = ~6iuj + uj')—lmg^^{uj; r, r'). (24) 



It is important to notice that deriving Eq. ()21|) and Eq. ()24|) with the help of Eq.( 
one has to assume that the temperature (entering Eq. (jl6j) ) is constant over the whole 
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space (where line 7^ 0). This is correct in thermal equihbrium, when all absorbing 
bodies are at the same temperature, and hence Eq. ()2H) provides the correlator for the 
equilibrium EM field. But it does not provide the correlator for the nonequilibrium 
thermal EM radiation, when there is radiation energy exchange between bodies with 
different temperatures. This correlator can be also calculated using Eq. (fTHjl . but the 
integral over the source point can not be simplified using Eq. lj^ . 



4. Dephasing by an infinite thick plate 

In what follows we consider a simple geometry when the dephasing body is a half- 
space 2; < with e{Y,u) = e{uj) and the two beams move in vacuum 2; > in a plane 
z = d parallel to the interface. In this geometry the Green function of the EM field is 
translational invariant in the x, y plane and hence can be presented as follows ^H] 

r, r') = 1 ^ e^'^(^-^')(?„;3(^, k\z, z'), (25) 

where R is the component of r in the {x, y) plane and k is a vector in this plane. Because 
of the special geometry we are interested in the Green function Q for z = z' = d and 
a,P = x,y, and will denote gai3{uj,'k\d,d) = gapitj,^^)- Using the explicit expressions 
for Qajs given in |16j, one can write 

gaf3iuj,k) = gtiuj.k) 

where 



fc2 2 



al3 



gi{uj,k), (26) 



9iA-^^ k) = 9i,t{^^ k)*, gai3iuj, -k) = gafsiuj, k) = gaf^i-^^, k)*. (27) 

The gai3{uj,]i) are related to the correlator of the tangential components of the electric 
field in the plane z = d. Using Eqs.(j2Zl) can check from Eq.ipTj) that for a, P = x,y 

{E4R)EAR'))^ = j d'k e^'^(^-^')(E„i?;3)^k (28) 

with 

(E^Ef^U = [-) coth— (-Im^7./3(a;,k)). (29) 

The classical current in beam 1 is ji(r, t) = e\i{t)6(R — Ili{t))6{z — d), where Ri(t) 
is the trajectory of beam 1 and vi(t) = dlii[t)/dt is the electron velocity in this beam. 
One can present 

ji(r, t) = e6{z -d)l^^ldu ™ (h)^, (30) 

where the radiation amplitude is 

(li)k.= / :^vi(t) e"^^W. (31) 
J 27r 

The relevant frequencies and wave vectors are those of the EM field created by the 
electron in beam 1. Similar expressions can be written for beam 2. We now introduce 
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into the dephasing integral, Eq.(jHI), the correlator T> expressed in terms of Q according 
to Eq. lj^ . and the Fourier expansions of ji 2 and Q according to Eq- ilHUj) and Eq.((2S)), 
we find, using Eqs.(j2ZI), 

2hc^ Jo 2^ 2f J (^^~^"^^"'3(w,k)} [{lahJhhc. + c.c], (32) 

where (l)ka; = (li)kaj — (l2)ka;- The contribution to this integral comes from frequencies 
and wave vectors which are present simultaneously in the fluctuation spectra g{u,\i) 
and the radiation spectra (l)kcj- Using Eq. (j^!)|) one can rewrite the integral over d'^k as 
follows 

1 f°° ( " ^ 

dkk -2lmgt{uj, A;)(|/k.r) - Mgiico, k) - gt{uj, A;)](|(l)k^H , (33) 



2vr JO 

where (...) means angular average and /k<^ = k(l)kaj/fc- 

For slow enough electrons one can use the dipole approximation (DA), when the 
term kRi(t) in Eq. (j31|) can be neglected. The condition for this is u; ^ A;f , where uj and 
k are the typical frequency and wave vector of the EM fluctuations contributing to the 
integral K, and v is the characteristic electron velocity. In this approximation (li)kaj = 
(li)^, and e(li)^ is the radiating dipole moment. Now (I)ka; = (li)a; — (12)0; = (1)(^- We 
substitute (1)^^ into Eq. (j32|l and as a result in the DA the dephasing integral is 

^ = 7T^ / — coth — S{uj) I (1)^|2, (34) 

where 



2nc2 /n 27r 2T 



|(1)J' = |(g.r + Kg.r, S{uj) = -Im / g,(uj,k). (35) 



d'^k 

S{lj) is related to average amplitude of the tangential electric field at distance d from 
the interface. From Eq.(j2HI) at R' = R one finds 



{EX = {Ei + Et^)^ = - [-) coth — S{uj). (36) 

It is convenient to represent K = Kp + Kg, where the two terms are the contributions 
to the integral of the domains, correspondingly, k < uj/c and k > uj/c. In the first 
domain the wave vector component perpendicular to the interface kz = a/ {oj/c)'^ — k'^ 
is real, which means that this domain corresponds to waves propagating perpendicular 
to the interface (PW), while in the second one kz is imaginary and it corresponds to 
evanescent waves (EW). 

5. Spectral densities 5f(c(j, /c) 

Using the exphcit expressions for gaii given in (TH], one can find (for > 0) 

g,,{u, k) = [e^P-(OF,,{0 - G,,(0], (37) 

where ^ = |k|/A;o with ko = uj/c, and 

^ ,1 ^ ,lv — u V — e(u;)u 

u u V + u V + e{uj)u 
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4;rcr CO 
Figure 1. Borderlines in the (lli, k) plane, see text. 

with u = [1 — i'^Y^'^, V = [e{u) — ^^]^/^,Im-u > 0, Imv > and p = 2kod. The upper 
sign corresponds to / and the lower to t. 

Consider the spectral densities {—Imgi^t) in the {uj, k) plane at > 0, A; > (see 
Figure^). One important borderline in this plane is = 1, i.e. k = uj/c, which, as noted 
above, separates the propagating waves (PW) domain below it from the evanescent 
waves (EW) domain above it. In the PW domain ^ < 1 and 

2tt 

- Img^ito, k) = —[Gi^tii) - ReFi,t(0 cosp^z(0], (39) 

fco 

while in the EW domain ^ > 1 and 

- lmgi^t{u:, k) = -^e-Pv^ReF,,,(0. (40) 

For e{uj) = 1 one finds u = v and F{^) = 0. This corresponds to empty space, in which 
case the spectral densities —Imgi^ti^^, k) ^ only in the PW domain, where in this case 

-lmgi,tiuj,k) = ^G^{0. (41) 

The same result is obviously obtained far from the interface, when c? — > oo, since one 
can neglect the oscillating or decaying term in Eq. ljHTj) . 

The second important borderline is = |e(a;)|^/^. For a non-dispersive lossless 
dielectric, Ime = 0, e = n"^ > 1, this borderline is simply k = uj/cn, where c„ = c/n 
is the light velocity in the dielectric. One easily finds that —Imgi^ti^^, k) = above the 
second borderline, for ^ > n. Below it, for 1 < ^ < n, one has 

ReF,,(0 = -^(e-a^/^ 

In the generic case of arbitrary complex e{uj) the spectral densities {—Imgi^t) are non-zero 
in the whole [u, k) plane. 

In what follows we will consider two cases: (i) a highly polarizable lossless dielectric, 
Ime = 0, e ^ 1, and (ii) a "good" conductor, Ime Ree ~ 1. In the last case we write 



± 1 



(42) 
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e{uj) = Co + i^^Tra/uj), where a is the conductivity, and assume that a is larger than all 
relevant frequencies. Both cases correspond to |e| ^ 1 and it is instructive therefore 
to investigate the limit |e| = oo. In this limit F{^) = G{C,) and the spectral densities 
—lmgi^t{uj, k) ^ only in the PW domain, where 

- Im^7^,,(^, k) = '^[l- cospu{0]Gi,t{0- (43) 

One can see that near the interface with an ideal conductor {a = oo) or an ideal 
dielectric (eo = oo) the tangential electric field fluctuations are suppressed. In the case 
of a conductor it is obvious from the boundary condition Et = 0. In the case of a 
dielectric this boundary condition is also effective, since Et would mean, due to the 
continuity of Et, infinite energy density in the dielectric or infinite displacement current. 

Now we turn to the spectral density corrections which are due to finite |e|. To 
investigate the role of these corrections we use the following expansions. Well below the 
second borderline, i.e. at ^ -C |e|^/^, one has 

FiiO = GiiO - 4e-i/2 + 0(6-1), ^ ^1 _ 2,-1] + 0{e-'/'). (44) 

Well above the second borderline, i.e. at ^ ^ kl^^^; one finds 

R,.(f)=>|^? + .(|^)^4 + 0(r'). (45) 

At ^ ~ |e|i/2 obviously F^{^) c^\e\^/^. 

It follows from Eq. ()44|) that in the PW domain the finite |e| corrections are small. 
However as we will see later, these corrections are important in the EW domain at small 
distances from the interface d. To account for these corrections in a dielectric one can 
use the explicit expressions given by Eq.(jl2I), but for a conductor the situation is more 
complicated. For a good conductor it is convenient to use the surface impedance C{uj) 
and the skin depth 6{u) defined as follows, 

C(^) = {u;/8nay/' ^ {2\e{u)\)-'/\ 6{u) = ^^J^y/, = 2C(^)^ • (46) 

In these terms the second borderline is k = \/2/5{lo) or ^ = {\/2C,{ijj))~^ . Since C,{uo) -C 1 
and 5{uj) <^ k^^ when uj <^ Ana, this borderline in the EW domain is well above the 
first borderline k = uj/c. 

In between the borderlines, oj / c <^ k <^ 5{ijj)~^, one can find using Eqs. pHl . 

ReF,(0 = -4C(^), ReFi(0 = -4C(^)'e (47) 

Above the upper borderline, k ^ 5{ijj)~^, one finds from Eq. (j45|) . using (e — l)/(e+ 1) = 
1 + 4i^(u;)^, the dominant term to be 

ReFz,i(0 = -K{^fi- (48) 

Now we find from Eq. ()40|) the spectral densities for a good conductor in the EW domain. 
In between the borderlines 

- lmgi{u, k) = ^CHe-''", -lmgt{uj, k) = ^Ci^fke-'^", (49) 
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while above the upper borderhne 

- Img^Au^, k) = ^aufke-"^', -lmg^,{u, k) = ^a^fke-'''^-^. (50) 

In the EW domain the fluctuations are small, because of the small surface impedance, 
and are strongly suppressed at > d~^, since random fields created by fluctuations of 
the random currents with wavelength shorter than d are averaged at distance d. 



6. Spectral density S{uj) 

The spectral density, which enters in the DA, can be split into contributions of the PW 
and EW domains, S{uj) = Sp{uj) + Se{uj), with 

S.iu) = - IfcoRe^'rfar^i^KO -GKO], (51) 

-| poo 

To avoid misunderstanding we note that the spectral density due to a half-space z < 
with e(a') calculated in Ref. [TTj is not for the equilibrium EM field, but is for the 
radiation into a zero temperature half space z > 0. 

In the limit |e| = cxd only PW contribute and one finds easily the spectral density 

S{u;) = Sp{uj) = - 
c 

At large distances from the interface, d ^ k^^, one has S{uj) = (2/3) /cq, which 
correspond to EM field fluctuations in empty space. (The factor 2/3 appear because 
only two tangential components of the electric field are relevant). At small distances, 
d <^ /cq ^, the fluctuations are suppressed, S{uj) ~ d^. For an ideal conductor or ideal 
dielectric the contribution to S{uj) comes from ^ — 1 (i.e. k c:^ u/c) independent of the 
value of the parameter p = 2kod. 

Now we turn to the finite |e| corrections. One can see from Eq. p^ that the 
corrections to Sp{u) are small, and one can use for Sp{uj) the result given by Eq.(jH2)). 
But this is not the case for Se{u)) when d is small and the decay of the exponent is 
this integral is slow. For a dielectric, Se{uj) can be calculated using Eq. ()42j) and when 
e = ^ 1 the result is 

2 11 

d<t: Xn- Se{uj) = -kon ; ci > A„ : Se{uj) = ^, (53) 

where A„ = Cn/uj = {kQn)~^ is the wave length in the dielectric. The first of Eqs.(jHSl) 
means simply that the fluctuations of Et at d A„ in the vacuum are the same as in the 
dielectric, which is consistent with the continuity of Et at the interface. Eqs.()53|) clearly 
demonstrate the importance of finite |e| corrections at small distances. Comparing Se{uj) 
for (i ^ A„ with Sp{uj) from Eq. (j3^ . one can see, that the ideal dielectric approach, 
when S{iu) is dominated by PW, is valid only at ci ^ Xn~^^^ = Xnn^^^, while at smaller 
distances S{u) is dominated by EW. According to the ideal dielectric approximation, 



2 
3 



cosp smp 



P 



(52) 
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Eq.(|^. near the interface the fluctuations are suppressed compared to free space, while 
from the first of Eqs. ()53|) it follows that they are enhanced compared to free space by 
a factor n. The contribution to the integral Se{uj) comes from k ~ A~^, when d <^ A„, 
and from k ~ when d ^ A„. When S'e(u;) dominates, both cases correspond to large 
imaginary k^. In other words, the fiuctuations contributing to Se{uj), are due to near 
fields localized close to the interface. 

The situation in a conductor is different, since from Eq. pHj) one can see that the 
integral Se{u!) diverges at d = 0. The near-field fiuctuations in the case of a conductor 
are 

d « 6{u;) : S.{u) = ^ = d» 6i.) : S^iuj) = ^ = (54) 

The near field/far field crossover point is d^ — {c/uj)C{ujy^^ and for all d the relevant 
k ~ d~^. When d <^ 5{uj) the relevant k are above the borderline k = \/2/5{uj) and 
the first of Eqs. ljK^ is obtained using Eq. ljHnjl . while when d ^ 6{u!) the relevant k are 
below this borderline and the second of Eqs.(|Mj) is obtained using Eq. ip^ . 

Comparing 5*6(0;) for a dielectric and a conductor one can see that they are similar 
not very close to the interface, if one replaces n^^ by C{uj) and A„ by S{uj). The behavior 
very close to the interface is different, since for a dielectric Se{uj) is finite, while for a 
conductor it diverges. (This singularity is cut-off if one takes into account the spatial 
dispersion of the conductor e pTj). 



7. Dephasing by near fields 

To simplify the picture of dephasing we assume that the electrons smoothly accelerate 
and decelerate. In other words, we assume that the electron motion has only one 
characteristic time scale r, which is the time of fiight from the emitter to the detector, 
and only one length scale L, which is the trajectory length. The characteristic electron 
velocity is defined as v = L/t. We also assume, following the experimental situation, 
that the beams are close, i.e. the distance between them a is small compared to L. The 
frequencies of the EM waves emitted by the electrons are of the order of u = r"^ and 
wavelengths are of the order of A = cr. In this model the angle averages entering the 
integrals in Eq. ljHHj) and Eq. ()H4j) can be presented as follows: 

where z = ujt/2, y = kL, and the small factor 6 = {a/LY appears because the beams 
are close and the effective current ji2 is smaller than the beam currents ji and j2. The 
functions \E' decay fast enough at > and k > L~^, restricting the integration in 
the {uj,k) plane in Eq. ljH^ within the rectangle □ = [ < u; < r~^, < A; < L~^], 
shown in Figure |21 by dashed lines. 

The frequency integral in Eq. ()34|1 contains three characteristic frequencies, namely, 
T/h (the frequency of the EM field thermal fiuctuations), (the frequency radiated by 
the electron), and c/d (the frequency which enters the spectral densities g). Assuming 



Decoherence of ele 



13 



k 




(ew) , k=a>/v 



L 



-1 



-1 



k = co/c^^ 
k = cojc 



-1 



-1 



CD 



T 



Figure 2. Integration domains in the (w, k) plane, see text. 



v/c = 0.1 and L = 10 cm, we have r = 3 x 10~^ s. The electro-optical system, 
which creates, guides and detects the beams, is at room temperature, and this is 
the temperature of the EM field surrounding the beams. At room temperature 
h/T = 2.5 X 10~^^ s and obviously always T/h ^ r~^. Hence one can replace in 
the integral Eq. ()32|) coth.{hLL!/2T) by its classical high temperature approximation. 

From the spectral densities calculated in secEland seclHlit follows that when |e| ^ 1 
the contribution Kp of the PW can be calculated as for an ideal conductor or dielectric. 
To calculate this contribution one can employ the DA, Eq. p^ . because when |e| = cxd 
the density S{uj) comes from k ~ uj/c, and for non-relativistic electrons the condition 
for the DA to be valid, namely kv <^ u, is satisfied. 

First we calculate the dephasing near an ideal conductor or ideal dielectric, given 
by Kp, Eq. ljMjl . Substituting there from Eq. (j33j) we obtain 



c?<A: Kp = hp2KQ—. 

Here Kq is the dephasing in free space, a = /he is the fine structure constant and the 
numerical factors 6 ~ 1 are given in terms of integrals 



The second of Eqs.(|3^ demonstrates that near an ideal conductor and dielectric the 
dephasing is suppressed. For the parameters used one obtains, neglecting numerical 
factors, Kq ~ 106'. This means that for well separated beams (0 1) at room 
temperature the dephasing due to thermal fiuctuations can be significant. In the existing 
experiments, however, the beams are very close, a = lOO/xm, and the dephasing is 
neghgible, Kq ~ 10 




(56) 




(57) 
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Now we turn to EW contribution K^, which is responsible for the enhancement of 
the dephasing near the interface. As one can see from Figure |21 this contribution is 
large for small electron velocities v, when the overlap of the rectangle □ with the EW 
domain is maximal. Motivated by this we consider first the simpler case of dielectric for 
f -C c„ = c/n. Since for a dielectric the spectral densities —lmgi^t{^,k) vanish above 
the borderline k = uj/cn, the integration domain overlaps only with the "bottom" of the 
rectangle □, where k is small, meaning that one can employ the DA. Being interested 
in an almost ideal dielectric (n ^ 1), we substitute the spectral density from Eqs. (|H!?|l 
into Eq. ljMj) and find the dephasing due to near fields to be 

< A„ : Ke = ^Kon; (58) 

3J 2 

> A„ : Ke = be Ko—r, be = -— -, 

nd^ 8 Jo 

where A„ = cr/n. Comparing the second of Eqs. ()58|1 with Eq. one finds the far 
field - near field crossover for dephasing to be dx — Xn~^^^, obviously the same as for 
S{u!). Contrary to the predictions of the ideal dielectric approximation (n = oo) the 
dephasing of beams moving near the interface is not suppressed compared to empty 
space, but enhanced by a large factor n/2. Since for most dielectrics n does not exceed 
10, the condition v <^ c/n is not very severe for non-relativistic electrons, but on the 
other hand, the enhancement of the dephasing near the interface is not very strong. 

The situation is much more complicated for conductors, since the spectral densities 
—Imgi^ti^, k) do not vanish above the upper borderline k = \/2/5{ijj) and the integration 
domain overlaps with the whole rectangle □. The parameter which plays the role of 
1/n in case of a conductor is C = (87r(Tr)^^/^, i.e. the surface impedance calculated for 
the characteristic frequency r^^. Copper at room temperature has a = 5 x 10^''sec~^ = 
6 X 10^ (ficm)"^, so for a good conductor this impedance can be as small as 10~^, and 
hence the restriction v <^ C,c can be very severe. As a result one has to consider velocities 
larger than (^c, when the DA might be invalid. Consequently the calculations are very 
involved, so we first present the results, discuss them, and sketch the calculations in 
the Appendix. In what follows we present the results for k = K/Kq and the crossover 
distance from near-field to far-field dephasing. The results are given in terms of the 
trajectory length L, the radiated wave length A = cr, the surface impedance (, and the 
skin depth 6 = 2X(. All numerical factors of order one are omitted. 

Three velocity intervals are relevant, namely, 

A: w/c<C; B: C < ^^/c < f^"; C : t'^<^v/c. (59) 

The crossover distances from near-field to far-field dephasing in these intervals are as 
follows 

A + B: dx=t'^X; C : d^ =t'^ {\^ / L). (60) 

The far-field dephasing, at d ^ dx, in all velocity intervals is given by Eq. ()5(ij) . The 
near-field dephasing, at d ^ dx, is different in different velocity intervals. 
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In interval A 

2 2 

d<^L: K = 6\/L^; L <^ d <^ 5 : K = 6X/d^; 5 < < (ix : K = SX/d^. (61) 
In interval B: 

d<^L: K = L<^d<^d^: k = 5\/d'^. (62) 

In interval C: 

d<^dy,: K = SX/L"^. (63) 

As one can see from the above results, in the velocity interval A + B the crossover dx 
is the same as for S{uj) and L <^ d^ X. In this velocity interval K depends on d 
in a non- monotonous way, reaching a minimum at d^, where /t ^ C^/^. In the velocity 
interval C one finds 5 -C fix L and approaching the interface K decays monotonously 
till d ~ (ix, where it saturates at k ~ ({v/c)~^. For all velocities K is finite at ci = 0, 
since at very small d the DA is invalid, and the singularity d~^ in S{u!) is cut-off by the 
ineffectiveness of wave vectors k > L~^. 

As was already mentioned, the dephasing Kq in empty space is very weak, and 
this is why the possible enhancement of K near the interface is of special interest. 
Looking for the ratio rj = K{d = 0)/Ko one can see from the above results that the 
dephasing near the interface is enhanced compared to that in empty space only for small 
enough electron velocities, when v/c <^ C^^^5 i-G- the interval A and in the smaller 
velocity part of interval B. Since in the experiment the fixed parameter is not r, but L, 
and hence ( depends on v, it is more convenient to use a different parameter, namely 
7 = c/SnaL. In terms of this parameter the velocity interval A is f/c <C 7 and the 
dephasing enhancement in this interval is 1] = 'j{v/c)~'^. The smaller velocity part of 
interval -B is 7 <^ f /c <^ 7 and here ^7 = 7 ^^'^{v/c)^'^^'^. (Note also, that the necessary 
condition C ^ 1 reduces to 7(f/c) ^ 1 and is always satisfied when 7^1.) For the 
parameters used above one finds 7 ~ 10^^° and 7 ~ 10^'^. It is clear now that in 
the case of a good metal the dephasing is enhanced only for relatively slow electrons 
and is not very high. For example, when v/c = 10~^ one finds r] ~ 10. Much stronger 
dephasing can be achieved with a high resistivity semiconductor, for example Si with 
(7 = 1 (f2cm)~^, in which case 7 ~ 10~^ and for v/c = 10~^ one finds rj ~ 10^. 

8. Dissipation versus dephasing 

The coherence of the electrons in the beams can be destroyed only if there are 
mechanisms which allow their energy to be dissipated. When there are no absorbing 
bodies in the EM environment of the beams, Eq.Q describes dephasing related to 
the dissipation of electron energy by radiation of EM waves "to infinity". In fact it 
means that the energy is dissipated in very far bodies, not included in the consideration 
explicitly. Eq.ijHJ is formally valid also when the beams are within a lossless cavity, if 
the small absorbtion in the walls is still large enough to prevent EM field buildup in the 
cavity. 
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If electrons in the two beams move along close trajectories and with similar 
velocities, ji2 = ji — j2 is small and the dephasing is weak. When the distance between 
the trajectories is small compared to the correlation length of the EM field in the 
direction perpendicular to the beams, the random electric fields in adjacent points of 
the two trajectories fluctuate synchronously, and as a result electrons in both beams 
change their phases also synchronously, which means that the beams remain mutually 
coherent. It does not mean, however, that the energy losses in the beams, defined by the 
currents ji and j2 separately, are small. There is one additional very important difference 
between dephasing and dissipation. Using the relation {D{aya''aD{a))T = n + 
where D{a) = exp[aa''' — a*a], one can prove that the energy radiated by a classical 
current into a thermal EM field is 

W = ^J dtdt'J drdrT{r,t)f{r',t') 2^K(r,t), A^(r',t')]- (64) 

Hence, in strong contrast to dephasing, the energy losses of the beam electrons do not 
depend on the environment temperature. 
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Appendix 

In what follows we sketch the calculations of the results presented in Eqs.lfUIH). 
and ()63|). There are two contributions to Ke, namely K'J, coming from above the 
borderline k = \/2/5{uj), and K'^, coming from between the borderlines k = uj/c and 
k = \/2/5{uj). These two contributions can be estimated using —lmgit{uj,k) from 
Eqs. fl5()|l and (^HI), correspondingly. One can also see from FigElthat when d ^ L the 
cut-off factor e~'^'"^ in —Im gi^t{u , k) selects from the rectangle □ only its "bottom" and 
hence the DA is valid. 

In the velocity interval A the lengths hierarchy is L ^ 5 ^ A. When d 5 the 
main contribution is K'J, where gi = gt, and the second term in Eq. ()33|) vanishes. As a 
result 

—2 

Ke = K = AiTo^ r % r dyy'e-^^'^'^^y^'riz, y). (65) 
Jo L-^ Jo Jo 

If d ^ L one can put z = 0, which corresponds to the DA, in agreement with what was 
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stated above, and obtain 



li d <^ L one can put d = and the integral is a numerical factor of order one. When 
d ^ 6 the main contribution is K'^ with gi dominating, and in addition the DA is valid. 
Substituting the second of Eqs. (j3^ into Eq.(|Mj) one finds 

In the velocity interval B + C the lengths hierarchy is 5 ^ L ^ A. Here the main 
contribution is always K'^ and one can check that gt can be neglected compared to gi. 
As a result 

If d ^ L one can put y = and obtain 

When d ^ L one can put c? = and the integral is a numerical factor. These results are 
valid in the whole velocity interval B + C . The separation appears when one compare 
the near-field and far-field contributions and finds that dx ^ L m B , while dy, L m 
C. 
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